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Dynamics  and  Control  of  a  Biomimetic  Vehicle  Using  Biased 
Wingbeat  Forcing  Functions:  Part  I  -  Aerodynamic  Model 
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An  aerodynamic  model,  for  a  minimally  actuated  flapping  wing  micro 
air  vehicle  (FWMAV),  is  derived  from  blade  element  theory.  The  vehicle 
considered  in  this  work  is  similar  to  the  Harvard  RoboFly,  except  that  it 
is  equipped  with  independently  actuated  wings.  A  blade  element-based 
approach  is  used  to  compute  both  instantaneous  and  cycle- averaged  forces 
and  moments  for  a  specific  type  of  wingbeat  motion  that  enables  nearly 
decoupled,  multi-degree-of- freedom  control  of  the  aircraft.  The  wing  po¬ 
sitions  are  controlled  using  oscillators  whose  frequencies  change  once  per 
wingbeat  cycle.  A  technique  is  introduced,  called  Split-Cycle  Constant- 
Period  Frequency  Modulation  with  Wing  Bias,  that  provides  a  high  level 
of  control  input  decoupling  for  vehicles  without  active  angle  of  attack  con¬ 
trol.  This  technique  allows  the  frequencies  of  the  upstroke  and  downstroke 
of  each  wing  to  differ  such  that  non-zero  cycle- averaged  drag  can  be  gen¬ 
erated.  Additionally,  a  wing  bias  term  has  been  added  to  the  wingbeat 
waveform  and  is  utilized  to  provide  pitching  moment  control.  With  this 
technique,  it  is  possible  to  achieve  five  degree-of-freedom  control  using  only 
two  physical  actuators.  The  present  paper  is  concerned  with  the  derivation 
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of  the  instantaneous  and  cycle- averaged  forces  and  moments  for  the  Split- 
Cycle  Constant-Period  Frequency  Modulation  with  Wing  Bias  technique. 
Implementation  of  the  wing  bias  is  discussed  and  modifications  to  the  wing- 
beat  forcing  function,  which  are  necessary  to  maintain  a  continuous  wing 
position,  are  made. 


I.  Introduction 

For  years,  researchers  have  been  interested  and  intrigued  by  the  flight  of  flapping  wing 
insects.  Some  of  the  earliest  research  efforts  on  flapping  wing  vehicle  flight  took  place  in 
the  1950’s  and  1960’s.1-4  Osborne1  used  blade-element  analysis  to  represent  the  forces  on 
a  flapping  wing.  Jensen2  studied  the  flight  of  locusts  while  Vogel3,4  analyzed  the  fruitfly 
Drosophila  and  calculated  a  lift  coefficient  for  this  insect.  More  recently,  flapping  wing  micro 
air  vehicles  (MAVs)  have  received  a  great  deal  of  interest  from  the  research  community  due 
to  their  potential  to  achieve  insect-like  maneuverability.5-11  The  ability  to  mimic  the  flight 
behavior  of  insects  could  enable  such  a  vehicle  to  perform  missions  which  larger,  fixed  wing 
vehicles  are  unable  to  perform,  such  as  intelligence,  surveillance,  and  reconnaissance  in  urban 
environments  and  indoor  locations.  These  potential  capabilities  have  prompted  continued 
research  into  flapping  wing  flight  dynamics,  control  law  development,12,13  and  fabrication.7,8 

Biologists  have  been  interested  in  insect  flight  for  some  time.  Ellington9, 14-18  provided 
detailed  analysis  of  hovering  insect  flight,  including  kinematics,  lift  and  power  requirements, 
and  aerodynamic  mechanisms.  Other  researchers  have  studied  the  flight  of  butterflies,19 
Diptera,20  bumblebees,10,11  and  the  hawkmoth  Manduca  Sexta.21-23  The  desire  was  to  char¬ 
acterize  and  understand  the  mechanisms  that  enable  insect  flight. 

The  aerodynamics  of  flapping  wing  flight  are  complex.9,14-18,24-2'  Wing-wing  interactions, 
wake  capture,  leading-edge  vortices,  wing  rotation  and  acceleration  are  examples  of  some  of 
the  complex  phenomena  that  occur.  Dickinson28  has  studied  the  unsteady  aerodynamics 
of  model  wings  at  low  Reynolds  numbers.  Sane  and  Dickinson6  have  performed  mineral  oil 
experiments  on  a  flapping  wing  and  generated  a  quasi-steady  aerodynamic  model  of  the  wing, 
while  Sitaraman29  developed  a  numerical  model  of  the  unsteady  aerodynamics  of  flapping 
flight. 

In  this  work,  a  simplified  aerodynamic  model,  which  does  not  take  into  account  unsteady 
effects,  is  utilized  for  control  law  design.  Deng,30,31  Doman, 12,32  and  Oppenheimer13  have 
developed  cycle-averaged  control  laws  for  flapping  wing  micro  air  vehicles.  These  control 
strategies  used  wing  flip  timing  and  mean  angle  of  attack  as  control  variables  or  included 
an  additional  control  effector  to  adjust  the  center-of-gravity  (CG)  of  the  vehicle.  The  con¬ 
trol  strategy  introduced  by  Deng30  required  4  actuators.  Oppenheimer13  utilized  split-cycle 
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constant  period  frequency  modulation  without  wing  bias  to  control  a  six  degree-of-freedom 
(DOF)  unsteady  simulation  model  of  a  FWMAV  using  three  actuators.  Two  actuators  con¬ 
trolled  wing  position  while  the  third  controlled  a  bobweight.  Here,  the  need  for  the  bobweight 
actuator  is  eliminated  by  modifying  the  wingbeat  profile  to  include  a  non-zero  bias.  The 
split-cycle  constant  period  frequency  technique12, 13  is  still  utilized,  however,  the  wingbeat 
can  be  biased,  which  has  the  effect  of  providing  control  of  the  pitching  moment.  Due  to 
the  periodic  nature  of  flapping  wing  motion,  the  forces  produced  by  the  wings  are  periodic 
functions  of  time.  Thus,  continuous  control,  of  the  forces  and  moments  produced  by  the 
wings,  is  not  an  option  and  control  variables  should  only  be  modified  at  the  start  of  a  cycle. 
Therefore,  a  cycle-averaged  control  law  is  developed  and  in  this  paper,  the  cycle-averaged 
forces  and  moments  for  a  flapping  wing  MAV  are  computed.  The  vehicle  used  in  this  work 
is  similar  to  the  RoboFly7  with  a  key  difference  being  that  the  present  vehicle  has  indepen¬ 
dently  actuated  wings.12  Wood7,8,33  has  shown  that  liftoff  can  be  achieved  with  these  types 
of  vehicles.  The  goal  of  this  work  is  to  provide  a  wingbeat  forcing  function  and  control  law 
that  can  yield  controlled  six  DOF  flight  of  the  vehicle  using  only  two  physical  actuators. 

A.  Preliminary  Assumptions 

The  purpose  of  this  investigation  is  the  development  of  a  control  strategy.  Hence,  an  idealized 
model  is  desired,  which  allows  one  to  study  the  dynamic  behavior  of  the  aircraft.  This 
objective  limits  the  fidelity  of  the  aerodynamic  and  structural  models  that  can  be  utilized. 
In  light  of  the  above  objective,  the  following  assumptions  are  made: 

•  There  are  no  aerodynamic  interactions  between  the  left  and  right  wings. 

•  There  are  no  aerodynamic  interactions  between  the  wings  and  the  fuselage. 

•  The  passive  wing  rotation  joint  is  on  a  limit  when  the  wing  angular  velocity  in  the 
stroke  plane  is  non-zero. 

•  The  2D  sectional  aerodynamic  coefficients  are  known  and  constant  throughout  each 
stroke.  Three  dimensional  and  unsteady  effects  from  leading  edge  vortex  and  wake 
phenomena  are  not  present. 

•  The  air  mass  surrounding  the  vehicle  is  quiescent. 

•  Aerodynamic  forces  and  moments  are  the  sole  result  of  wing  motion,  i.e.,  only  flight 
conditions  in  the  neighborhood  of  hover  are  considered. 

•  The  bandwidth  of  the  piezoelectric  actuators  exceeds  the  wingbeat  frequency. 
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A  reduced-order  model  for  control  law  design  is  developed.  Ultimatley,  the  controller 
developed  will  be  tested  on  a  truth  model  that  relaxes  some  of  the  assumptions  made  above 
and  includes  unsteady  aerodynamic  effects.  Although  the  unsteady  model  is  not  developed 
here,  it  will  be  a  focus  of  future  research  efforts. 

II.  Split-Cycle  Constant-Period  Frequency  Modulation  With 

Wing  Bias 

The  wingbeat  position  command  structure  presented  by  Doman12  was  shown  to  be  ca¬ 
pable  of  producing  horizontal  and  vertical  forces  as  well  as  rolling  and  yawing  moments. 
The  inability  of  that  control  strategy  to  independently  generate  pitching  moments  led  to  the 
addition  of  a  bobweight  actuator.  Here,  a  wing  position  bias  is  introduced  that  enables  the 
cycle-averaged  center-of-pressure  of  the  beating  wings  to  be  moved  such  that  cycle-averaged 
pitching  moments  can  be  manipulated  without  a  bobweight.  It  will  be  shown  that  the  new 
control  strategy  is  capable  of  generating  six  independent  forces  and  moments  using  only  two 
physical  actuators. 


4> uLW{t )  —  ALw  cos  [(ulw  —  SLW )  t]  +  r}LW 
4>dLW{t)  =  Alw  cos  [( uLW  +  &lw)  t  +  £lw]  +  Vlw 


and 


< Krw  (t)  ~  Arw  cos  {(ivrw  ~  &rw)  t]  +  rjRW 
<t>dRW{t)  =  Arw  COS  [(urw  +  <JRw)t  +  Crw]  +  V RW 

where  (pULW(t)  is  the  wing  position  of  the  left  wing  on  the  upstroke,  <f>dLW(t )  is  the  wing 
position  of  the  left  wing  on  the  downstroke,  c t>URW  (t)  is  the  wing  position  of  the  right  wing  on 
the  upstroke,  4>dRW(t)  is  the  wing  position  of  the  right  wing  on  the  downstroke,  t]rw,  rjLw  are 
the  wing  bias  terms,  and  ARW,  ALW  are  the  wingbeat  amplitudes.  Practical  implementation 
requires  that  the  bias  terms  be  delayed  by  one  full  wingbeat  cycle  for  reasons  that  will 
become  evident  shortly.  Changing  the  amplitude  will  have  a  similar  effect  as  changing  the 
fundamental  wingbeat  frequency  because  these  variables  change  the  wing  velocities  in  the 
same  way.  Hence,  either  the  amplitude  or  the  wingbeat  frequency  can  be  a  control  variable. 
In  this  case,  left  and  right  wingbeat  frequency  will  be  used  for  control  purposes,  while  the 
amplitude  is  maintained  constant  (with  the  exception  of  the  necessary  last  quarter  stroke 
amplitude  change  for  wing  position  continuity  with  nonzero  bias,  which  will  be  discussed 
shortly). 

Figure  1  shows  the  effects,  on  the  wingbeat  forcing  function,  of  changing  the  wing  bias. 
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The  wing  bias  is  forced  to  be  constant  over  each  full  cycle.  Since  rj0  ^  rji  ^  772,  discontinuities 


0  0.5  1  1.5  2  2.5  3 

Normalized  Time 


Figure  1.  Wing  Position  With  Variable  77. 

exist  in  the  wing  position.  This  is  a  situation  that  cannot  physically  exist.  To  overcome  this 
issue,  the  second  half  of  the  downstroke  is  modified  by  the  addition  of  a  cosine  wave.  For 
the  first  half  of  the  downstroke,  the  wing  position  forcing  function  is  exactly  as  shown  in 
Equations  1  and  2.  For  the  second  half  of  the  downstroke,  the  wing  position  forcing  functions 
are  modified  to  the  following  form 

<t>dLW{t)  =  Alw  cos  [{ujlw  +  aLw)  t  +  £lw]  +  Vlw  +  AALw  cos  [( u>lw  +  &lw)  t  +  Clw]  (3) 
and 

4>dRw  (t)  =  Arw  cos  [{urw  +  &rw)  t  +  £rw]  +  Vrw  +  ^ Arw  cos  [(ujrw  +  aRW)  t  +  £Rw]  (4) 

The  terms  A  Alw  and  A  ARw  are  used  to  maintain  a  continuous  wingbeat  position  forcing 
function.  Figure  2  shows  the  results  of  adjusting  the  last  half  cycle  of  the  downstroke.  In  this 
case,  the  wing  position  is  continuous.  Because  frequency  was  chosen  as  a  control  variable 
instead  of  amplitude,  this  work  utilizes  constant  amplitudes,  i.e.,  Alw  =  ARw ■  Allowing 
the  amplitude  to  vary  between  the  up  and  downstrokes  would  force  changes  to  the  wingbeat 
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functions  at  the  half  and  full  cycle  points  of  the  wing  motion  to  maintain  a  continuous 
wing  position,  whereas,  changing  the  wing  bias  at  the  end  of  a  full  cycle  only  requires  one 
modification  to  maintain  a  continuous  wing  position. 


Wing  Position  vs.  Time 


Figure  2.  Wing  Position  With  Variable  77  Altered  to  Maintain  Continuity. 

Figure  3  shows  the  points  in  the  wingbeat  cycle  where  04,  Si,  r/i,  A A*  are  calculated  and 
applied.  Note  that  these  variables  are  calculated  at  the  beginning  of  the  upstroke  and  that 
AAj  is  applied  for  the  last  half  of  the  downstroke  to  ensure  wing  position  continuity.  The  bias 
term  is  delayed  by  one  full  wingbeat  cycle  and  then  is  applied  for  the  next  complete  wingbeat 
cycle.  Thus,  at  the  beginning  of  an  upstroke,  u 4  and  h*  are  calculated  and  immediately  applied 
for  a  complete  cycle.  Also,  AAj,  Aj,  and  r/*  are  calculated  with  A  A*  being  applied  for  the 
last  half  of  the  current  wingbeat  downstroke  and  A,  and  r/j  being  applied  for  the  next  full 
cycle. 

I11  order  to  calculate  the  amplitude  adjustment  terms,  consider  the  end  of  a  complete 
wingbeat  and,  for  brevity,  consider  only  the  left  wing  (right  wing  calculations  are  similar). 
Using  Figure  3  and  Equation  3,  it  can  be  seen  that  at  the  end  of  the  Erst  complete  cycle, 
when  t  —  — ,  the  modified  downstroke  wing  position  is 

4>dLW(t)\t=2^  =  ALWq  +  r)LWo  +  AALWl  (5) 
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Wing  Position  vs.  Time 


Figure  3.  Points  In  Wingbeat  Cycle  Where  Control  Parameters  Are  Calculated  And  Applied. 

while  at  the  beginning  of  the  second  cycle  upstroke,  the  wing  position,  Equation  1,  is  given 
by 

4>uLW{t)\t=^.  =  ALWi  +  VLWi  (6) 

Equating  the  expression  in  Equations  5  and  6  and  solving  for  A ALWl  yields 

AALWl  =  ALWl  +  r]LWl  —  ( ALWq  +  r]LWo)  (V 

In  general,  this  expression  becomes 

AAj  =  Ai  +  rji  —  {Ai_i  +  (8) 

Thus,  in  order  to  determine  AAt,  the  system  needs  information  about  the  amplitude  and 
bias  during  the  next  wingbeat  cycle. 
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The  forcing  functions  driving  the  positions  of  the  wings  are 


<t>uLW(t )  =  Alw  cos  [(ijJlw  -  Slw)  t]  +  TjLW  (upstroke) 


,(t)  =  Alw  cos  [(ulw  +  &lw)  t  +  C.LW }  +  Vlw  ^first  i  of  downstroke^ 

,(t)  =  ( Alw  +  A Alw)  cos  [(u )lw  +  &lw)  t  +  £lw\  +  Vlw  f  second  ^  of  downstroke 


4>uRW(t)  =  Arw  cos  [{urw  -  Srw)  t]  +  r/RW  (upstroke) 

<PdRW  (t)  =  Arw  cos  [(urw  +  cjrw)  t  +  £rw]  +  r]RW  ^first  ^  of  downstroke 


<f>dRW(t)  =  ( Arw  +  A  Arw)  cos  [{urw  +  aRW)  t  +  Cmv]  +  Vrw  (second  -  of  downstroke 


ft  is  important  to  point  out  that,  in  Equations  9  and  10,  the  time  index  subscripts  on  the 
variables  have  not  been  included.  It  should  be  understood  that  for  wingbeat  cycle  i,  the 
biases  and  amplitudes  used  are  r]Lwi-1^rlRWi-1i  ARWi  l,  Arw^  while  the  fundamental  fre¬ 
quency  urwi^lw,  split-cycle  parameters  Srw,  Slw,  and  amplitude  modification  parameters 
AArw,AAlw  are  the  current  time  values.  These  variables  are  produced  by  the  cycle- 
averaged  control  law  presented  in  Part  II.34  In  other  words,  bias  and  amplitude  are  delayed 
by  one  complete  wingbeat  cycle  prior  to  their  use  in  generating  the  wingbeat  forcing  func¬ 
tions. 


III.  Instantaneous  Blade  Element  Model  of  Flapping  Wing 
Vehicle  with  Independently  Actuated  Wings 

In  this  section,  an  aerodynamic  model  of  the  flapping  wing  vehicle  is  developed.  One  of 
the  objectives  of  this  work  is  to  develop  a  model  which  yields  insight  into  the  controllability 
of  this  class  of  vehicles.  In  light  of  this  objective,  the  aerodynamic  model  is  void  of  any 
unsteady  effects  and  is  based  solely  on  blade  element  theory.  The  use  of  computational 
fluid  dynamics  or  finite  element  models  would  provide  a  higher  fidelity  model,  however, 
these  techniques  would  not  allow  designers  to  make  vehicle  modifications  and  rapidly  asses 
their  effects  and  have  a  tendency  to  mask  the  most  pertinent  factors  in  force  and  moment 
generation. 
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A.  Instantaneous  Aerodynamic  Forces  and  Centers-of-Pressure  in  Wing  Plan- 
form  Frames 

The  aerodynamic  forces  are  derived,  using  blade  element  theory,  for  triangular  shaped  wings 
that  have  two  degrees  of  freedom,  namely,  angular  displacement,  </>(£),  about  the  wing  root 
in  the  stroke  plane,  and  angular  displacement  of  the  planform  about  the  passive  rotation 
hinge  joint,  which  is  equivalent  to  wing  angle-of-attack,  a,  in  still  air.  The  lift  and  drag  due 
to  the  wings  was  previously  calculated12  and  is  summarized  here. 

The  lift  and  drag  can  be  expressed  as  the  product  of  time  invariant  parameters  and  time 
varying  functions 


L  =  kL(j){t )2 
D  =  kDc/)(t)2 


(11) 


where 


kL  =  P~ CL(a)IA 
kD  =  P- CD(a)IA 


p  is  the  atmospheric  density,  I  a  is  the  area  moment  of  inertia  of  the  planform  about  the 
wing  root  and  Cl  (a),  Cd(cx)  are  the  lift  and  drag  coefficients,  which  are  computed  using  the 
following  expression  from  Sane  and  Dickinson5 


CL(a)  =  1.58  sin  (2.13a  -  7.2)  +  0.225 
CD(a )  =  -1.55  cos  (2.04a  -  9.82)  +  1.92 


(13) 


The  only  variable  that  can  be  actively  manipulated  to  control  the  instantaneous  aerodynamic 
forces  is  the  angular  velocity  of  the  wing,  (j){t).  Differentiating  Equations  9  and  10  yields 


<j>ULW(t)  =  -Alw  ( uLW  -  SLW )  sin  [{uLW  -  SLW)  t]  (upstroke) 

<t>dLW{t)  =  ~Alw  (ujlw  +  aLW )  sin  [( uLW  +  aLW)  t  +  £ LW }  ^first  ^  of  downstroke^j 
4>dLW(t)  =  —  ( Alw  +  A Alw)  (lvlw  +  Tlvu) 


(14) 


*  sm 


[(wlw  +  ctlw)  t  +  £lw\  ^second  -  of  downstrokej 
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and 


&uRW(t)  —  -ARW  (urw  —  $rw )  sin  [{urw  —  Srw)  t]  (upstroke) 

4>dRW  (t)  =  -Arw  (urw  +  vrw)  sin  [(ujrW  +  aRW)  t  +  £RW]  ( first  ]-  of  downstroke  J 

V  J  (15) 

YdRW(t)  —  ~  (Arw  +  AArW)  (u>rw  +  &rw) 

*  sin  [( ujrw  +  (Jrw)  t  +  £rw}  ^second  -  of  downstroke^ 

The  instantaneous  lift  and  drag  of  the  left  and  right  wings  on  the  upstrokes  and  downstrokes 
are  obtained  by  substituting  Equations  14  and  15  into  Equation  11  to  yield 

LUlw  =  kLA2LW  (u:LW  -  5lw)2  sin2  [(uLW  -  SLW )  t]  (upstroke) 

LdLW  =  kLA2LW  (uLW  +  aLW)2  sin2  [(uLW  +  aLW )  t  +  £LW]  ^first  ^  of  downstroke^) 

LdLW  =  kL  (Alw  +  A Alw)2  (ulw  +  ctlwY 


*  sin2  [( ulw  +  &lw)  t  +  Clw]  ^second  -  of  downstrokej 
Dulw  =  kDA2LW  (uLW  -  5lw)2  sin2  [(uLW  -  hw)  t]  (upstroke) 

DdLW  =  kDA2LW  ( uLW  +  aLW)2  sin2  [( uLW  +  crLW )  t  +  £ LW }  ^first  ^  of  downstroke^) 
DdLW  =  kD  ( Alw  +  A Alw)2  ( ujlw  +  &lw)2 

*  sin2  [(ujlw  +  &lw)  t  +  Clw]  ^second  -  of  downstroke^ 


(16) 


(17) 


Lurw  =  kLA2RW  (uRW  -  5rw)  sin2  [(ujrw  -  SRW )  t]  (upstroke) 

LdRW  =  kLA2RW  (uRW  +  aRW)2  sin2  [(ojrw  +  aRW)  t  +  £RW]  ^first  ^  of  downstroke^) 

LdRW  =  kL  ( ArW  +  AArwY  (ujrw  +  &rwY 

*  sin2  [( ujrw  +  &rw)  t  +  £,rw]  ^second  -  of  downstroke^ 

Durw  =  kDA2RW  (uRW  -  SrW )2  sin2  [(uRW  -  SRW )  t]  (upstroke) 

DdRW  =  kDA2RW  (ujrw  +  ctrw)2  sin2  [(uRW  +  aRW)  t  +  £RW]  ^first  i  of  downstroke^) 

DdRW  =  kn  (Arw  +  AARW)2  (wrw  +  &rwY 

*  sin2  [(ujrw  +  cr rw)  t  +  Crw]  ^second  -  of  downstroke^ 

With  the  relationships  between  the  body,  roots,  spars,  upstroke  planform  and  downstroke 


Approved  for  public  release;  distribution  unlimited. 
10  of  34 


planform  axis  systems  established  by  Doman,  et.  al., 12  the  instantaneous  values  of  lift  and 
drag  on  each  wing  may  be  transformed  into  the  body-axis  coordinate  frame.  If  the  air  mass 
is  quiescent,  then  the  relative  wind  is  parallel  to  the  stroke  plane,  which  is  coincident  with 
the  x-y  plane  of  the  local  spar  frames.  Therefore,  the  lift  and  drag  forces  are  conveniently 
expressed  in  the  spar  coordinate  frames,  which  can  be  transformed  to  the  body  frame  using 
rotation  matrices.12  The  aerodynamic  forces,  associated  with  each  wing  and  stroke  expressed 
in  the  body  frame,  are  summarized  in  Table  1. 


Force 

Body  Frame 

RW  Upstroke 

Fs  = 

URW 

Lurw 

Durw  sRl  (^) 

.  DUrw  COS  (t) 

RW  Downstroke 

F  f  = 

dRW 

LdRW 

$dRw  (^) 

_  DdRW  COS  (t)dRW(f) . 

LW  Upstroke 

pB  — 

ULW 

ULW 

DulW  &ULW  (f) 
PuLW  CO®  (fouLW  (^)_ 

LW  Downstroke 

■p-B 
r  d.LW 

LdLW 

-DdLW  sin  4>dLW  it) 
—  DdLW  cos  <pdLW(t)_ 

Table  1.  Aerodynamic  forces  expressed  in  local  spar  and  body  frames. 

In  Table  1,  the  subscripts  denote  which  wing  and  stroke  while  the  superscript  defines  the 
axis  system  in  which  the  force  vector  is  written.  For  example,  the  subscript  Urw  denotes 
the  right  wing  upstroke  while  d^w  denotes  the  right  wing  downstroke.  The  superscript  B 
denotes  the  body  frame. 

The  center-of-pressure  of  each  wing  is  conveniently  expressed  in  the  local  wing  planform 
frame  associated  with  each  stroke.  The  centers-of-pressure12  of  each  wing  on  the  upstroke 
and  downstroke  are  summarized  in  Table  2.  In  Table  2,  Ar^  is  the  position  vector  from  the 
center-of-gravity  to  the  origin  of  the  right  wing  coordinate  system,  Arf  is  the  position  vector 
from  the  center-of-gravity  to  the  origin  of  the  left  wing  coordinate  system,  and  x^pp)y^pp 
define  the  location  of  the  center-of-pressure  in  the  local  wing  planform  frame.  These  position 
vectors  are  explicitly  expressed  as 
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CP  Location 

Body  Frame  Expression 

RW  Upstroke 

r  B  = 

CP  Urw 

x^pp  sin  a  +  Axp 

x™p  sin  (f)URW  cos  a  +  y™p  cos  (j)URW  +  f 
-a^P  cos  (pURW  cos  a  +  yWp  sin  +  Az% 

RW  Downstroke 

B  _ 

IcP  dRW 

x^pp  sin  a  +  Aar| 

-x™p  sin  (t)dRW  cos  a  +  y ™p  cos  +  f 

_x^p  cos  <; bdRW  cos  a  +  ycHp  p  sin  <f>dRW  +  Azp_ 

LW  Upstroke 

r  B  = 

cPuLW 

x^p  sin  a  +  Axf 

-Xcp  sill  (pULW  cos  a  -  y™  cos  (pULW  -  ^ 
-x^p  cos  cos  a  +  t/^p  sin  (j)ULW  +  Az£  _ 

LW  Downstroke 

r  B  - 
1  CP  dLW 

a;^p  sin  a  +  Axf 

%cp  sm  <PdLW  COS  a  -  y™  cos  (j)dLW  -  j 
_x^p  cos  cos  a  +  t/Jp  sin  (j)dLW  +  Azf  _ 

Table  2.  Centers-of-pressure  expressed  in  body  frame. 


Arf  = 


A  xx 


w_ 

2 


Arf  = 


Axf  -f 


Azf 


(18) 


where  w  is  the  width  of  the  vehicle. 


B.  Aerodynamic  Moments  in  Body  Frame 

The  expressions  for  the  aerodynamic  moments  associated  with  each  wing  and  stroke  are 
given  by12 


=  rc  p  x  Ff 

URW  cP  Urw  Urw 

lvrB  —  r  B  v 

d-RW  cPdRW  ditw 
ML  =  rcp uLW  X  <W 

^-dLW  =  rcPdiw  X 


(19) 
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Carrying  out  the  cross  product  operations  and  substituting  the  values  from  Tables  2  and  1 
into  the  expressions  in  Equation  19  yields12 

^URW  \VcP  ~b  "2  COS  0URW  "b  AzR 
urw  \Vcp  ^urw  "b  A~^>]  Du rwAxR  COS  4>urw  •  •  • 

[LUrw  cos  a  +  DUrw  sin  a]  x^p  cos  (j>URW  }  (20) 

{  ^URW^XR  ^URW  —  LURW  \_~2  "b  Vcp  ^URw]  •  •  • 

[LUrw  cos  a  +  DUrw  sin  a]  x sin  (bURW  } 

D<1rw  \llcp  ~b  ~2  cos  fidRw  "b  a%r  sin 
JdRw  \_Vcp  tfidRw  "b  A~^>]  “b  DdRWAx r  COS  (pdRW  T  •  •  • 

Mfw  =  |  [Aw  COS  a  +  sin  a]  x™p  cos  (pdrtw  }  (21) 

{DdRWAx ft  sin  4>dRW  L dRW  \_~2  "b  Vcp  cos<pdRW\  +  . . . 

[LdRW  cos  a  +  DdRW  sin  a]  x^p  sin  (j)dRW  } 

Aw  [~yZP  ~  f  COS(PuLW  ~  Az£  sin  <j>ULW] 

{ Lulw  \VcpP  sin  fiiiLw  +  Azl]  -  DulwaXl  cos  (j)ULW  -  ... 

[LUlw  cos  a  +  DUlw  sin  a]  x^p  cos  <t>ULW  }  (22) 

{DUlwAxl  sin  4*ulw  "b  LUlw  ["2  “b  HCp  COS  (fiuRwl  +  ■■■ 

[LUlw  cos  a  +  D*lw  Sln  a]  X^F  Sln  } 

KP  +  f  cos  (j)dLW  +  A z£  sin  0diW] 

{  KP  sin  ^dLW  +  A  z£]  +  DdLWAxf  cos  (j)dLW  +  ... 

MdLW  =  [LdLW  cos  a  +  sin  a]  x, ™p  cos  }  (23) 

{-AwAxf  sin  4>dLW  +  LdiW  [f  +  t/^P  cos  -  . . . 

l  [  cos  a  +  Ddw  sin  a]  x)Jp  sin  (f>dLW  } 

Equations  20-23  provide  the  expressions  for  the  instantaneous  aerodynamic  moments  gener¬ 
ated  by  each  wing  at  any  point  in  a  wingbeat  cycle  in  the  body-axis  coordinate  frame. 

In  this  work,  the  relationships  between  the  cycle-averaged  forces  and  moments  and  the 
wing  kinematic  parameters  that  can  be  manipulated  for  control  are  derived.  These  pa¬ 
rameters  are  culw,  tu rw ,  8lw?  $rw,  Vlw ,  and  Vrw,  which  are  the  fundamental  frequency, 
split-cycle  parameters,  and  wing  biases  for  the  left  and  right  wings,  respectively,  are  derived. 
Feedback  control  laws  based  on  cycle-averaged  forces  and  moments  will  allow  a  vehicle  to 
track  desired  angular  and  spatial  positions  in  a  mean  sense;  however,  because  of  the  true 
periodic  nature  of  the  aerodynamic  forces,  the  vehicle  will  exhibit  limit  cycle  behavior  in  a 
neighborhood  about  the  mean  position. 
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IV.  Definite  Integrals  for  Computation  of  Cycle- Averaged 
Aerodynamic  Forces  and  Moments 

In  this  section,  expressions  for  the  compute  cycle-averaged  forces  and  moments  are  de¬ 
rived.  Since  the  introduction  of  the  split-cycle  parameter  changes  the  frequency  of  the  cosine 
wave  when  (u j  —  5)t  =  7r  in  each  cycle  and  because  the  downstroke  is  split  into  two  parts,  it 
is  convenient  to  split  the  integrals,  used  to  evaluate  the  cycle-averaged  forces  and  moments, 
as  follows 

G{<j>d(t))dt  +  /  G(<j>d(t))dt  (24) 

J  7r(3iu  — 2d) 

2cj  (cu  —  S) 

where  the  downstroke  is  split  into  first-half  and  second-half  parts.  Note  that  G  in  Equation  24 
is  a  generalized  force  that  represents  either  a  force  or  moment. 

In  order  to  calculate  the  cycle-averaged  forces  and  moments,  it  will  be  necessary  to 
evaluate  numerous  integrals.  Many  have  no  indefinite  integral  solutions.  For  example,  many 
of  the  integrands  will  be  of  the  form  cos(coscuf)  or  sin  (cos  cut).  Fortunately,  definite  integrals 
involving  such  functions  exist  over  the  intervals  of  interest  for  the  present  problem  and  can 
be  derived  from  results  presented  in  Gradshteyn  &  Ryzhik.35  The  solution  of  many  of  these 
definite  integrals  involve  a  Bessel  function  of  the  first  kind,  Ji(-),  as  well  as  a  Struve  function, 
H i(-).35  For  convenience,  the  solutions  to  the  definite  integrals  that  are  required  to  compute 
the  cycle-averaged  forces  and  moments  are  provided. 

For  the  upstroke,  there  are  three  integrals  of  interest.  The  first  is 

7T 

h  =  f  sin2  [(cu  -  S)t]  dt  =  -  (25) 

Jo  2  (a;  o) 

The  second  integral  used  during  the  upstroke  is 

7T 

/2  =  f  sin2  [(a;  —  5)t]  sin  (A  cos  [(a;  —  S)t]  +  rj)  dt  (26) 

Jo 

Using  the  following  identity, 

sin  (  fi  +  rj)  =  sin  f3  cos  rj  +  cos  (3  sin  rj  (27) 


2?r 


G(4>u(t))dt  + 


^  7r(3cu  — 25) 
2ct>(u j  —  5) 
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Equation  26  can  be  rewritten  as 


/2  =  /  sin2  [(tu  —  5)t]  sin  ( A  cos  [(cu  —  <5)t])  cos  ifdt 

Jo 

7 r 

f  u-S 

+  /  sin"  [(a;  —  5)t\  cos  {A  cos  [(a;  —  5)t ])  sin  rfdt  =  I21  +  I22 


Evaluating  the  integrals  yields 


hi  =  0  I2 2  =  ^Ji(A)  sin?7 


where  Ji(A)  is  a  bessel  function  of  the  Erst  kind.  The  last  upstroke  integral  necessary  to 
compute  cycle-averaged  forces  and  moments  is 

7T 

I3  =  /  sin2  [(a;  —  8)t]  cos  (A  cos  [(a;  —  S)t]  +  rj)  dt.  (30) 


Using  the  following  identity, 


cos  (/ 3  +  rj)  =  cos  (3  cos  rj  —  sin  f3  sin  rj 


Equation  30  can  be  rewritten  as 


^  /  u>  —  5 

I3—  sin"  [(a;  —  S)t]  cos  (A cos  [(a>  —  8)t])  cosrfdt 
Jo 

/E 

uj  —  8 

sin"  [(a>  —  S)t]  sin  [A  cos  [(a;  —  8)t])  sin  rfdt  =  J3 1  +  J32 


Evaluating  the  integrals  yields 

^31  =  A(w-5)  J^A)  cos  if  I32  =  0  (33) 

Since  the  second  half  of  the  downstroke  is  modified  to  ensure  wing  position  continuity, 
integrals  over  the  downstroke  become  more  complicated  than  those  over  the  upstroke.  In 
fact,  the  integrals  are  split  into  two  parts.  The  first  downstroke  integral  is 


7t(3cl>  —  28)  2-7T 

f  2oj(co-S)  f  — 

h=  sin2  [(a;  +  a)t  +  £]  dt  +  /  sin2  [(a;  +  a)t  +  £]  dt 

I  7T  /  7r(3ct)  —  26) 

LO-S  2lo(lu-6) 

_  71 -  7T  _ 

—  T 7 - 7  H — 77 - 7  —  -*41  +  1 42 

4  (u  +  &)  4  (u  +  cr) 
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The  fifth  integral  is 


7t(3cj  —  25) 

A  f  2lu(uj-8) 

h  =  sin2  [(w  +  o)t  +  £]  sin  (A  cos  [(a;  +  a)t  +  £]  +  rj)  dt 


+  /  sin"  [(a;  +  a)t  +  £]  sin  ({A  +  AA}  cos  [(a;  +  a)t  +  £]  +  rj)  dt 

I  7r(3a;  —  25) 

J  2lu(uj-8) 


Using  the  identity  in  Equation  27,  J5  can  be  written  as 


7r(3o>  —  2(5) 

A  f  2u> (uj  —  5) 

h  —  sin2  [(w  +  cr)t  +  £]  sin  ( A  cos  [(a;  +  a)t  +  £])  cos  rjdt 

uj  —  8 

7r(3a;  —  25) 

/2uj{uj  —  8) 

sin"  [(a;  +  a)t  +  £]  cos  ( A  cos  [{uj  +  a)t  +  £])  sin  rjdt 

2tt 

+  /  sin2  [(a;  +  a)t  +  £]  sin  ({A  +  AA}  cos  [(a;  +  a)t  +  £])  cos  rjdt 

2uj(uj  —  8) 

+  /  sin2  [{uj  +  a)t  +  £]  cos  ({A  +  AA}  cos  [(a;  +  a)t  +  £])  sin  7]dt 

J  7r(3a;  — 2 8) 

2oj(u)—8) 

=  hi  +  h2  +  -^53  +  hi 

where 

J51  =  2A^T)HM)  COS  71 
J52  = 

^53  =  2{A+AA}(u+<j)Hi{A  +  A  A)  COST? 

-^54  =  2{A+AA}(w+CT)  S^n  d 

and  ifi (•)  is  a  Struve  function.  The  last  integral  of  interest  is 


tt(3uj  —  28) 
a  I*  2uj(uj  —  8) 

h=  sin2  [(a;  +  cr)t  +  £]  cos  ( A  cos  [(a;  +  cr)t  +  £]  +  rj)  dt 


uj  —  8 
„  2ir 


+  /  sin2  [(a;  +  a)t  +  £]  cos  ({A  +  AA}  cos  [(a;  +  a)t  +  £]  +  rj)  dt 

I  7T (3cu  —  28) 

2uj(uj  —  8) 
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Using  the  identity  in  Equation  31,  Equation  38  can  be  rewritten  as 


7r(3it->  — 2(5) 

^  I*  2u>(lo  —  8) 

J6  =  /  sin"  [(a;  +  cr)t  +  £]  cos  (A  cos  [(a;  +  a)t  +  £])  cos  rjdt. 


^  —  28) 
2a }(u  —  8) 


o-8 
„  2tt 


sin2  [(a;  +  cr)t  +  £]  sin  (A  cos  [(a;  +  <j)t  +  £])  sin  rjdt 


+  /  sin"  [(a>  +  <r)t  +  £]  cos  ({A  +  AA}  cos  [(a;  +  a)t  +  £])  cos  rjdt 

I  7t(3cj  —  25) 

2a>(a>  — (5) 

„  2t r 


f  7r(3a;-2(5) 
2a;(cj  — (5) 


sin2  [(a;  +  a)t  +  £]  sin  ({A  +  AA}  cos  [(a;  +  a)t  +  £])  sin  rjdt 


—  hi  +  h2  +  -^63  +  hi 


(39) 


where 


J61  =  2A^Ta)JM)  COS?? 
hi  =  2A^Ta)HM)^d 

h3  =  2{A+AA}{uj+a )  h{A  +  A  A)  COS  7] 

hi  —  2{A+AA}(uj+a )  H\{A  +  A  A)  sin  77 


(40) 


V.  Cycle- Averaged  Forces  and  Moments 

The  cycle-averaged  aerodynamic  forces  and  moments  in  the  body  frame  will  now  be 
calculated.  The  cycle-averaged  forces  and  moments  will  eventually  form  the  basis  for  the 
computation  of  sensitivities  to  changes  in  control  parameters.  These  cycle-averaged  sensitiv¬ 
ities  form  the  control  effectiveness  matrix,  which  is  used  to  compute  the  wingbeat  kinematic 
parameters.  These  kinematic  parameters  produce  the  cycle-averaged  forces  and  moments 
that  are  commanded  by  a  feedback  control  law  that  regulates  the  mean  motion  of  the  vehi¬ 
cle. 


A.  X  Force 


Substituting  the  expression  for  the  instantaneous  x-body  force  from  Table  1  into  Equation  24 
produces 


—B  _  UJrw 
A  RW  ~ 


2t r 


wRW~SRW 


URW 


(■ t)dt 


.  n(3UJRW~2SRw) 
2ulRW  d’RW  ~  SRW  ) 


JdRW 


(' t)dt 


dRW  —°RW 


4"  /  /,  M  LdRW(t) 

I  n(3ulRW-2dRW ) 

2ujrw(u,rw~srw) 

Approved  for  public  release;  distribution  unlimited. 
17  of  34 


(41) 


Note  that  the  limits  of  integration  in  Equation  41  for  the  last  two  integrals  define  which  part 
of  the  downstroke  is  under  consideration  and  hence,  which  (pdIiW  (t)  from  Equation  10  to  use. 

The  integral  from  - —  to  n^3uRw~2SRw'>  js  the  first  half  cycle  of  the  downstroke,  while 

the  integral  from  n(3uJRw-2SRw)  |-0  _2zl.  jg  ^jie  seconq  }ia]f'  Qf  the  downstroke.  Substituting 

0  ^rw(Urw~  ORW)  Urw  0 

Equation  17  into  Equation  41  yields 


^ RW 


k'LURW 
2i r 


t3RW~°RW 


arw(uRW  —  &rw)  sin  \(^rw  ~  ^Rw)t\dt 


7r{3ulRW~2SRw) 

+  [  A2Rw(,URW  +  aRw)2  Avl2[(uRW  +  (?Rw)t  +  ^Rw\dt 


tlRW~°RW 

2  7T 


f  n(3uiIiW-25RW) 
2ulRW  (" RW  ~  5RW  ) 


(Arw  +  AARW)  (ujrw  +  vrw)  sin  [(ujrw  +  <?Rw)t  +  ^Rw]dt 


(42) 


Noting  that  the  time  varying  functions  under  the  integral  signs  are  of  the  form  I\  and  J4  as 
given  by  Equations  25  and  34  respectively,  Equation  42  can  be  written  as 


4  =  [A  rw(uRW  ~  Srw)2Ii  +  A2rw(ujrW  +  ctrwYIai 


2ir 


+kL  ( Arw  +  A  Arw)2  (ujrw  +  &rw)2I±2 


or  simply 


— b  _  k^RW 

A  RW  —  A 


«2  [ .  .  r  n,  ,  (^RW  +  <7rw) 


Arw  ( URW  —  Srw)  + 


{Arw  +  (Arw  +  AArw)2'\ 


Following  a  similar  procedure  for  the  left  wing,  it  can  be  shown  that 

A\w  (ulw  ~  $lw)  +  -  —  0 — {A2lw  +  (Alw  +  A Alw )2} 


— b  _  kLUJLW 

A  LW  ~  a 


(43) 


(44) 


(45) 


_ Q  _ jTf 

Note  that  both  XRW  and  X LW  are  positive  quantities.  At  a  hover  condition,  where  the 
x-body  axis  is  normal  to  the  surface  of  the  earth,  the  forces  produced  by  both  wings  act 
to  counter  the  weight  of  the  vehicle.  Also,  from  Equations  44  and  45,  it  can  be  seen  that 
wing  bias  has  a  secondary  effect  on  the  vertical  force  produced  and  only  appears  in  the 
AArw,AAlw  terms.  In  practice,  it  is  expected  that  A ARW  -C  1  and  A ALW  -C  1  and  thus 
minor  changes  to  vertical  force  will  be  experienced  by  a  non-zero  wing  bias.  Recall,  from 
Doman,12  that  when  no  wing  bias  is  used,  XRW  =  kLUJ^w  (2u>rw  —  Srw  +  ctrw)  and  similarly 
for  the  left  wing.  Comparing  this  previous  result  with  the  expression  in  Equation  44  and 
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setting  Arw  =  1  rad  shows  that  the  increase  in  x-body  axis  force,  that  results  from  the 
introduction  of  a  wing  bias,  is  given  by  kLU}RW(ujRW+t7RW'>  {2 A ARW  +  A A2RW}  and  similarly 
for  the  left  wing. 

B.  Y  Force 

The  cycle-averaged  y-body  force  for  each  wing  is  now  considered.  Substituting  the  expression 
for  the  instantaneous  y-body  aerodynamic  force  from  Table  1  into  Equation  24  produces 


^-pB  _  UJrw 
1  RW  ~ 


2vr 


U1RW^SRW 


n(3uiRW-2SRW) 

/  2uirw{^rw-5rw)  „  ,  n 

/  ^  UdRW  (t) 

sin[^w(t)]dt  + 

r- 

r 

1  ^{3uRW~2SRw) 

DdRW(t)  sin [<j>dRW(t)]dt 

“ RW~SRW 

2u1RW  (" RW  -  SRW  ) 

(46) 


Substituting  Equations  10  and  17  into  Equation  46  yields 

—A2Rw(uFtw  —  Srw)2  sin2[(o;Rvv  —  iw)t] 


—R  kDURW  f“RW-sRW  a2  /  x  \  2  •  2  r 


^  RW  ~ 


2vr  J0 

*  sin  (ARW  cos[(curvi/  —  Sr w)t\  +  i]Rw)  dt 

n  (3uirw-2Srw) 

+  °2 I  Rn(UJRW  Rw)  A2RW(uRW  +  ( trw )2  sin2 [{ujrw  +  °RW)t  +  £RW] 


dRW  ~°RW 


*  sin  (Arw  cos[(curvu  +  <?Rw)t  +  £rw]  +  Vrw)  dt 

,  kDURW  ,  A  A  ^ 


2tT  J  n(3u,RW~25Rw) 

2wRW  (“ RW~SRW  ) 


(Arw  +  A Arw)  ( u>Rw  +  (Trw)  sin  [(uRw  +  <?Rw)t  +  £rw] 


*  sin  ({Arw  +  A ARW}  cos[(u  rw  +  (?Rw)t  +  £rw]  +  Vrw)  dt 

Using  Equations  28,  29,  35,  36,  and  37,  Equation  47  becomes 
7b  koojRw 


(47) 


Y  rw  ~ 


2tt 


Arw  (urw  ~  Srw)  (hi  +  I22 )  +  Arw  (ujrw  +  <?rw)  (^51  +  h2) 


+  {Arw  +  AArw}2  (ujrw  +  &rwY  (I53  +  ^54) 


(48) 
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Simplifying  Equation  48  yields 

— I^dAr-w Ji{Arw)urw  0 ,jJRW  —  3rw)  sinrjRw 


yb  — 

1  rw  ~ 


ko^RW  {^>RW  +  °~Rw) 
4 

kD^RW  i^RW  +  CTrw) 


{Arw  (Ji(Arw)  sin Vrw  ~  Hi{Arw)  cos? Irw)} 


(■ Arw  +  AARW) 


■h  {Arw  +  AARW )  sin  r)RW 


(49) 


+  Hi(Arw  +  A Arw)  cos  riRW 


From  Table  1,  it  can  be  seen  that  the  left  wing  side  force  has  the  same  form  of  the  right 
wing  side  force,  except  for  a  minus  sign.  Hence,  for  the  left  wing,  it  can  be  shown  that 


^tB  _  koALW  Ji(Alw)u>lw  (a>lw  ~  &lw)  sin  r)LW 

*  LW  ~  o 


krjijJLW  (y>lw  +  o~lw) 
4 

kn^LW  {^lw  +  ®lw) 


{Alw  (Ji(Alw)  Anr]LW  —  Hi{Alw)  cos rjLW)} 


(■ Alw  +  As.Alw) 


J\{Alw  +  AAlw)  sinr)LW 


(50) 


+  Hi(Alw  +  A  Alw)  cos  t/lw 


Note  that  when  the  wing  bias  is  zero  [rjRw  =  0,  r/LW  =  0)  the  amplitude  continuity  term  is 

_ q  _ q 

also  zero  (A ARW  =  0,  A ALW  =  0)  and  Y RW \mw=o  =  0,  Y lw\Vlw=0  =  0.  Without  wing  bias, 
a  side  force  could  not  be  generated  even  when  using  split-cycle  frequency  modulation.  Now, 
the  wing  bias  has  added  the  ability  to  generate  side  forces,  although  for  small  wing  biases, 
it  is  expected  that  only  small  side  forces  could  be  produced. 


C.  Z  Force 

The  calculation  of  the  cycle-averaged  force  in  the  z-body  direction  follows  a  similar  procedure. 
Substituting  the  expression  for  the  instantaneous  z-body  aerodynamic  force  from  Table  1  into 
Equation  24  gives 


~yB  _  UJRW 
LRW  — 


2tt 


" RW~SRW 


DURW(t)  cos [<t>URW{t)]dt 


.  n{iulRW^2SRw) 
2ulRW  (“flW  ~SRW  ) 


-DdRW  it)  COS [<j>dHW  (t)\dt  + 


VRW~°RW 


J  7r(3ojjiW -2Srw) 

2u1RW  (“RW  -  SRW  ) 


-DdRW(t)cos[<pdRW(t)]dt 

(51) 
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Substituting  Equations  10  and  17  into  Equation  51  yields 


-yB  _  koOJRW 
Zrw  ~ 


xJRW~°RW 


A2Rw(uRW  —  SRW)2  sin2[(o;Rvu  —  dRW)t] 


*  cos  ( Arw  cos[(ujRw  —  SRw)t]  +  i]Rw)  dt 

/ni3ulRW-2SRw) 

2ulRW  (" RW  ~  SRW )  -2  /  \  2  •  2lY  \  ,  j.  “1 

ARW\URW  +  <7Rw)  Slir[(ufijy  +  <7 RW)t  +  ? RW\ 


‘1RW~°RW 


*  cos  (Arw  cos[(toflo'  +  Vrw)1  +  One]  +  ’/me)  dt 


(52) 


I  tv(2>uj -25 RW) 

2uirw{^rw~srw) 


(Arw  +  AArw)  (uRw  +  ctrw)  sin  [(uRw  +  <7  Rw)t  +  €rw] 


*  cos  ({Arw  +  AARW}  cos[(cur^  +  <7Rw)t  +  £ Rw ]  +  Vrw )  dt. 


Noting  that  the  terms  of  the  integrand  are  of  the  form  J3  and  J6  as  given  by  Equations  30 
and  38  respectively,  Equation  52  can  be  written  as 


_  k.DOJRW 

Zrw~^T 


A2Rw(uRW  ~  $rw)2  (hi  +  ^32)  —  A2RW(uRw  +  7Rw)2  (hi  +  dm) 


—  {ARw  +  AArw}2  (uRW  +  <7Rw)2  (hz  +  hi) 

Simplifying  Equation  53  yields 

koARw  Ji(ARw)uirw  (urw  —  dRw)  cos  r/Rw 


yB  — 
^  RW  — 


knuJRw  (urw  +  < 7rw ) 
4 

kD^RW  (OJRW  +  7Rw) 


{Arw  (Ji(ARw)  cos  r]Rw  +  Hi(Arw)  sin r)Rw)} 


(ARw  +  AARW) 


Ji(ARw  +  AARW)  cos  r/nw 


—  H\(ARw  +  AARw)  sin rjRw 


(53) 


(54) 
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For  the  left  wing,  the  cycle-averaged  z  force  becomes 


-p^b  _  koALW  Ji(Alw)ljlw  (ujlw  ~  Slw)  cos  r]LW 

Z  LW  ~  o 


k'D^LW  LW  +  &Lw) 

4 

ko^LW  LW  +  ffw) 

4 


{^Lir  (Ji(Alw)  cosrjLW  +  H1(ALW)  sin  rjLW)} 


{Alw  +  As.Alw) 


Jl{ALW  +  AAlw)  COS  1JLW 


(55) 


—  Hi(Alw  +  AALW)  sin r/LW 


The  result  for  the  z-body  axis  is  important  because  it  means  that  non-zero  cycle-averaged 
forces  can  be  generated  and  used  to  induce  fore  and  aft  linear  accelerations  by  symmetrically 
varying  the  split-cycle  parameters.  It  can  also  be  used  to  generate  rolling  moments  by 
asymmetrically  (right  to  left)  varying  the  split-cycle  parameters. 

The  results  for  the  cycle-averaged  forces  can  be  simplified  to  those  obtained  in  Doman12 
by  setting  rjRW  =  ijLW  =  0  and  Anw  =  ALW  =  1. 


D.  Rolling  Moment 

For  the  rolling  moment,  the  expression  in  Equation  24  becomes 


wB  -  URW 

1V1XRW  ~ 


VRW~°RW 


2vr 


'  tt(3ujrW-26rw) 
2u1RW  i^RW  ~SRW ) 


MXURW(t)dt  + 
MxdRW(t)dt  I 


.  k^RW-^Rw) 
2wrw{ujrw-srw) 


MxdRW{t)dt 


dRW  —°RW 


(56) 
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Substituting  Equations  10  and  17  into  Equations  20  and  21,  Equation  56  becomes 


B  _  OJ rw  ,  WP 


Mt.  rw  - 


koVcp  (uRW  —  6rw)  A 


2a2 
^ RW 


UJRW~SRW 


sin2[(oduu  —  $Rw)t]dt 


+  kp>—(ujRw  —  dRw)2  A2rw 

*  sim[(cjRVK  —  dRw)t]dt 


cos  { Arw  cos  [( oJRw  —  SRw)t]  +  Vrw} 


+  kpAzRA2RW(ujRW  —  SRw)' 


sin  {Arw  cos  [{(jJrw  —  dRw)t\  +  Vrw} 


*  sin2[(cuRn/  —  SRw)t]dt 


—  kD(uJRW  +  aRW )  ycp  Arw 


jt(3ujrw-2Srw) 


‘>RW\u1RW-6RW)  .  2 


sin  \(ujRw  +  &Rw)t  +  £rw]  dt 


—  koi^Rw  +  aRw)2-^A2RW 


' w(3ujRW~2SRW ) 


urw^rw~°rw)  .  2 


sin  [{wRw  +  crRW)t  +  imv] 


*  cos  {Arw  cos  [{urw  +  crRw)t  +  £rw]  +  Vrw}  dt 


ko{uRw  +  crRw)  AzrArw 


A3  “RW~25RW) 

2  a  B  a2  I  2uRw(uRW~sRw)  .  2 


sin  [(ujrw  +  &Rw)t  +  £rw] 


*  sin  {Arw  cos  [{ijJrw  +  Rw)t  +  £rw\  +  Vrw}  dt. 


kp(^Rw  +  ORw)2yZP  {Arw  +  AARW 


o  /  “ RW 

)  \ 

I  n(3uRW~2SRw) 
2ujrw(u1rw-srw  , 


sin2  [( ujRw  +  &Rw)t  +  £,rw]  dt 


~  kD(uRW  +  ctrw)277  {Arw  +  AArw ) 


2  /  “RW 


RW  ~25 Rw) 


° RW  \“ RW  RW  ) 


sin2  [( ijJRw  +  &Rw)t  +  £ rw ] 


*  cos  { (yARw  +  A Arw)  cos  [( ojRw  +  &R,w)t  +  £rw]  +  Vrw}  dt 


—  kp(uRw  +  ctrw)2Azr  {Arw  +  AARW) 


2  /  URW 

J  n(3ulRW-2SRw) 
2ulRW  (.uRW  -  SRW 


sin2  [{wrw  +  crRW)t  +  £rw] 


*  sin  {{Arw  +  AArw)  cos  [{urw  +  &Rw)t  +  £,rw]  +  Vrw}  dt. 
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Note  that  the  integrals  are  of  the  form  given  in  Equations  25  -  40;  thus,  Equation  57  can  be 


written  as 


mxrw  -  D  { (urw  -  Srw)2A2rw  y™1  /,  +  A zRI2  +  ~h  } 

- D  ^  aRw)2A2RW  y^pPI 41  +  A zR  (J51  +  152)  +  —  (hi  +  I&2)  } 

- ~^TT  D  {(URW  aRwy  (Arw  +  AArW)2  y^PI 42  +  Azr  (I53  +  154)  +  —  (^63  +  -^4) 


Substituting  the  results  for  the  definite  integrals  and  simplifying  yields 

— — B  kDARW  Urw  ( Urw  -  Srw )  r  WP  .  ,  .  w  .  0  A  B  ■  11 

Mxrw  — - ^ -  [Vcp  aRW  +  Ji(Arw)  {w  cos r]RW  +  2 A zR  sin  (jRw}\ 

kpARWURW  (urw  +  °~Rw)  WP  A  ,  17/4  \  ,  tt  (  a  \  •  1 

- - -  t/cp  +  w \J\(ArW)  cosijrw  +  Hi(Arw)  smrjRw} 


+  2A zR  {J\(Arw)  sin rjRW  —  Hi(Arw)  cos t/rw} 


kp  (Arw  +  AArw)  ujrw  (urw  +  <7rw)  wp  r  a  ,  a  a  \ 
- g -  Vcp  \arw  +  ^arw) 

+  w  {  J\  (Arw  +  A  Arw)  cos  qRW  —  Hi  (ARW  +  A  ARW)  sin  r)RW} 


kp  (Arw  +  A  Arw)  ojrw  (&rw  +  &rw) 


Hi  (Arw  +  A  Arw)  cos  i]rw} 


2  AzR  {  J\  (Arw  +  A  Arw)  sin  t/rw 
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Following  a  similar  procedure  for  the  left  wing,  it  can  be  shown  that 


Mxlw  ~ 


kDALWuLW  (uLW  -  8lw)  r  WP  A  i  t  t  A  \  f  ,  0  a  b  ■  n 

- - -  [ycp  Alw  +  Ji(Alw)  (re  cos  r]LW  +  2A zL  sm  rjLW\\ 


^dAlW^LW  \UlW  +  clw) 


yZPALW  +  w  {Ji(Alw)  cos i]LW  +  Hi(Alw )  smrjLW} 


2 («/i(^4lw)  shi i]lw  ~  Hi(Alw )  cos t/lw} 
ko  (Alw  +  A Alw)  u>lw  {^lw  +  &lw) 


ycp  {Alw  +  A  Alw) 


+  iv  {  J\  ( Alw  +  AAlw)  cos  t)lw  —  H\  ( Alw  +  AAlw)  sin  tjlw} 
.  ko  {Alw  +  A Alw)  ujlw  ( uLW  +  (Jlw) 


2A zL  {Ji  {Alw  +  A Alw)  sin  t]LW 


Hi  (Alw  +  A Alw)  cos  r]LW} 


(60) 


When  the  wing  bias  terms  were  not  included,  then  split-cycle  frequency  modulation  was 
necessary  to  generate  non-zero  cycle-averaged  rolling  moments.  Note  that,  in  this  case, 
without  split-cycle  frequency  modulation,  i.e.,  5RW  =  aRW  =  5LW  =  aLW  =  0,  it  would  be 
possible  to  generate  non-zero  cycle-averaged  rolling  moments  on  this  aircraft  with  wing  bias. 
Hence,  wing  bias  has  effectively  added  an  additional  two  control  parameters  (left  and  right 
wing  bias)  capable  of  generating  rolling  moments. 


E.  Pitching  Moment 

The  cycle-averaged  aerodynamic  moments  in  the  y-body  frame  (pitching  moment)  is  now 
calculated.  For  the  pitching  moment,  the  expression  in  Equation  24  becomes 


M, 


VRW 


VRW 

2vr 


" RW~SRW 


MVuRW(t)dt  + 


tt(3uirw-2Srw) 
2u1RW  {wRW  -  SRW  ) 


MydRW^)dt 


°RW~°RW 


+ 


/  ,  .  My. 

'  tt(  3ujrw-25rw)  yaRW 

2ulRW  (“ RW  -  5RW  ) 


(t)dt 


(61) 
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S3  tO 


Let  C\  =  (corw  ~  $Rw)t  and  C2  =  (uj rw  +  &rw)  t  +  £rw ■  Then,  substituting  Equations  10  and  17 
into  Equation  61,  while  using  Equations  20  and  21  yields 


w 


URW  A2 
2? r 


Arw(uRW  ~  Srw)  <  kpAzj) 


B  /  UJRW~°RW  .  2 


sin  C\dt 


—  kpAx^  /  Rn  sin2  C\  cos  (Arw  cos  C\  +  tjrw)  dt 


7  wp  ■ 

—  KDXCp  sm  a 


7  WP 

—  klxcp  cos  a 


sin2  C\  cos  (A cos  C\  +  tjrw)  dt 

_ 7T _ 

JRW~SRW  .  2  /a  ^  \  1 1 

sm  Ci  cos  (Arw  cos  Gi  +  t]rw )  dt 


+kLy , 


WP  I  UJRW~aRW  .  2 

cp 


sin2  Ci  sin  ( Arw  cos  Ci  +  tjrw)  dt 


+  -^-A2rw(ujrw  +  (JRW )2  {  kLAz% 


B  /  2^rwCrw~srw)  ■  2  ^  i± 

'  sm  C2  dt 


^RW~SRW 


T  kp>AxR 


n{3ulRW-2SRw) 

B  I  2u,RwCrW~sRw)  ■  2 


sin2  C2  cos  ( Arw  cos  C2  +  t]rw)  dt 


" RW~SRW 

tt(3uirw-2SrW) 


7  WP 

+  klxcp  cos  a 


7  wp  ■ 

+  kdxcp  sm  a 


2uirwCrw-5rw)  ■  2  ri  1  a  n  ,  \  n 

sin  C2  COS  {Arw  cos  C2  +  TjRw)  dt 


llRW~0RW 

n(3ujRW-26RW) 

2^rwCrw~srw)  .  2  ri  7  A  n  ,  \  n 

Sill  C'2  COS  (  Arw  cos  C'2  +  ?Jrw)  dt 


+  kLy , 


UJRW~°RW 
’r(3ulRW-2SRw) 

WP  /  2uRwipRW~5Rw)  o.^2 

cp  / 


sin2  C2  sin  (Arw  cos  C2  +  r/®v)  dt 


lJRW~°RW 


+  ™  (Arw  +  A  Arw)2  (wrw  +  orwY  {  kpAzp 


'  ■k(3ujFiW-28rw) 
2u1RW  CrW  -  SRW  ) 


sin2  C2dt 


,  »  b  /  “w 

T  kpAxR 


^{,3u,rw-25rw) 
2ujrw  Crw~srw  ) 

_  2?r 


sin2  C2  cos  ((Arw  +  AA^^)  cos  C2  +  tjrw)  dt 


7  vkp 

+  klxcv  cos  a 


cp 


7  WP  ■ 

+  kdxcp  sm  a 


I  ^(3ujrw  ~2SRw) 
2uirwCrw~5rw) 

r— 

I  " RW 

I  tt(3i*jrw-26rw ) 
2ulRW  CrW  -  SRW  ) 


sin2  C2  cos  ((Tuiy  +  A  Arw)  cos  C2  +  tjrw)  dt 


sin2  C2  cos  ((Arw  +  A. Arw)  cos  C2  +  tjrw )  dt 


+  kLy , 


vy  p  /  "w 

CP  /  7r(3i^flXy-2i5flvlA) 

2ujrwCrw~6rw) 


sin2  C2  sin  ((d^  +  A  Arw)  cos  C2  +  ijrw)  dt 


(62) 
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Note  that  the  integrals  are  of  the  form  given  in  Equations  25  -  40;  thus,  Equation  62  can  be  written 


My  aw  -  Srw)2  kLAzRh  -  kD  AxRI3  -  kLx™p  cos  al3 

-  kDx^p  sin  al3  +  kLy^,pI2  +  A2rw(ujrw  +  ctrw)2  kLAzRhi 
+  kf)Axp  (Iqi  +  h2)  +  kL3w>P  cos  a  (hi  +  I&2) 

+  kDx^pp  sin  a  (hi  +  1(2)  +  kLy^pp  (hi  +  I52)  (63) 

+  (Arw  +  A^Arw)2  (urw  +  vrw)2  kLAzPh2 
+  kr>  Axp  (1(3  +  -^64)  +  kLx%p  cos  a  (1(3  +  -^64) 

+  kDx^pp  sin  a  (I63  +  hi)  +  kLy^p  (/53  +  J54)  | 

Substituting  the  results  for  the  definite  integrals,  I\  -  hi  and  simplifying  produces 

■jtb  -Arw^RW  (urw  ~  $Rw)  /  wPi  ,  r  •  WP  .  a  1  \  T  ( A  \ 

Mvrw  = - 5 -  (cos  axcp  kL  +  |smaxcp  +  AxR)  kD)  J\(ARW)  cos  r)RW 

wp 7  7  /  A  \  ■  A-zRkLARw 

~  Vcp  kL-J\  (Arw)  sm  riRW - ^ - 

H - 4Arw -  C°S  aXcp  kL  +  ^  S1R  aXcp  +  AXr  *  kf>  ArW  Jl  ^Arw  ' cos  VRW 

+  Hi(Arw)  sin  tjrw  }  +  Arw  (Arw  +  A  Arw) 

*  {Ji  ( Arw  +  A  Arw)  cos  t/rw  —  Hi  (Arw  +  A  Arw)  sin  tjrw} 

.  WPi  URW  (URW  +  ctrw)  r,2  IT  (A  ^  •  TT  !  A  \  I 

+  ycp  kL - 4  Arw -  Arw  Wi(Arw)  sm  tjrw  -  H\(ArW)  cos  tjrw  j 

+Arw  (Arw  +  A  Arw)  {J\  (Arw  +  A  Arw)  sin  tjrw  +  H\  (Arw  +  A  Arw)  cos  i/rw} 

+  MilZHK pEL ±  IMl  [a%w  +  (Arw  +  A ARwf 

(64) 
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Following  a  similar  procedure  for  the  left  wing,  it  can  be  shown  that 


— — B  -ALWu>LW  {uLW  ~  hw)  /  WPi  ,  r  •  H'P  .  A  ,  \  t  <  a  \ 

My lW  = - ^ -  (c°s  axcp  kL  +  { sin  axcp  +  AxL  }  kD)  (ALW )  cos  r/nr 

WPi  t  :  a  \  •  Az^kiALw 

~  Vcp  kLJi(ALW)  sm 7ILW - - - 

,  ULW  (^LW  +  &Lw)  (  WPr  ,  f  •  WP  I  A  Bill  4  2  [  I  l  A  ^ 

H - - (cos  axcp  +  {  sin  axcp  +  AxLj£:B)  Alw  {Ji{Alw)  cosrlw 

+  H\{Alw)  smr)LW}  +  ALW  [ALW  +  A ALW) 

*  { Ji  ( Alw  +  A>.Alw)  cos  tjlw  -  Hi  {Alw  +  A^Alw)  sin  t/lw} 

.  WP i  ULW  {ULW  +  &LW)  [,2  rr/A  \  •  it  t  a  \  1 

+  ycp  - 4  Alw -  ALW  {Ji(ALW)suir]Lw  ~  Hi(Alw)  cos  tjlw} 

+Alw  ( Alw  +  AAlw)  {Ji  {Alw  +  AAlw)  sin  t/lw  +  Hi  [Alw  +  J^Alw)  cos  tjlw} 

+  ^^LW^LW  +  ^LW)  ,Aiw  +  {Alw  + 

(65) 


Note  that  without  split-cycle  frequency  modulation,  i.e.,  Srw  =  grw  =  $lw  =  aLW  =  0,  there 
still  exists  a  non-zero  cycle- averaged  pitching  moment  if  A zp  %  0  and  A zp  %  0.  Setting  A z^  = 
0,  A zp  =  0  will  yield  a  zero  cycle- averaged  pitching  moment  when  the  split-cycle  parameters  and 
wing  biases  are  zero,  which  is  a  desirable  feature  for  maintaining  hover.  Also  note  that  the  wing 
bias  appears  throughout  these  expressions  and  thus  will  have  an  impact  on  the  pitching  moment  of 
the  vehicle.  To  obtain  a  more  quantitative  estimate  of  cycle- averaged  pitching  moment  change  due 
to  wing  bias,  Equation  64  can  be  rewritten  with  A z^  =  0,  sin  rrw  =  0,  cos  tjrw  =  1,  Arw  =  1  rad, 
and  [cos  ax^ppkL  +  {sin  ax^pp  +  Ax{j}  ko)  =  M  to  obtain 


-B  —urw  {urw  ~  &rw) 


MJi(l) 


Ww  {urw  +  (?rw) 


M  [Ji(l)  +  (1  +  A Arw)  Ji  (1  +  A^Arw)] 


WPr,  urw  [lorw  +  vrw) 


+  yZ  kL 


[—Hi[l)  +  (1  +  A Arw)  H i  (1  +  A Arw)\ 


When  AA  =  0,  Ti(l)  +  (1  +  A Abw)  J\  (1  +  AAbw)  =  2«7i(l)  =  .8801  and 

— i?i(l)+(l  +  A  Arw)  H\  (1  +  A  Arw)  =  0.  For  AA  =  0.2  rad ,  Ji(l)+(1  +  A  Arw)  Ji  (1  +  A  Arw)  = 
Ji(l)  +  1.2Ji(1.2)  =  1.038  and  —H\{1)  +  (1  +  0.2)  H\  (1.2)  =  0.1344.  Thus,  a  change  in  the  second 
and  third  terms  in  Equation  66  on  the  order  of  15%  can  be  expected;  therefore,  even  a  small  wing 
bias  can  significantly  alter  the  cycle- averaged  pitching  moment.  Figure  4  shows  Hi  (A)  and  Ji(A) 
vs.  A.  The  slopes  of  the  Bessel  function  of  the  first  kind  and  the  Struve  function  about  A  =  1  are 
nearly  the  same,  thus  it  is  expected  that  both  the  second  and  third  terms  in  Equation  66  will  im¬ 
pact  the  cycle-averaged  pitching  moment.  In  other  words,  when  computing  the  change  in  pitching 
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moment  due  to  a  change  in  ?/.  the  chain  rule  is  used  to  obtain  =  c)JiIqaa^W'> 

and  dHlly]^l^RW^>  =  dHl^AA^W^  9/driL™  •  Since  the  derivatives  of  the  Bessel  function  of  the  first 
kind  and  the  Struve  function  with  respect  to  A  are  non-zero  about  A  =  1,  the  second  and  third 


terms  in  Equation  66  will  create  pitching  moment  increments. 


Struve  Function,  H^A),  and  Bessel  function  of  first  kind,  J^A),  \s.  A 


Figure  4.  H\(A)  and  J\(A)  vs.  A. 


F.  Yawing  Moment 

The  cycle- averaged  yawing  moment  (z-body  moment)  is  calculated  in  the  same  fashion.  Substitut¬ 
ing  the  expression  for  the  instantaneous  z-body  moment  from  Equations  20  and  21  into  Equation  24 
produces 


M- 


-B 


ZRW 


URW 

2tt 


URW~°RW 


M. 


ZURW 


(■ t)dt  + 


/. 


t(3u. 


VRW{UJRW~°RW  ) 


M. 


zdh 


,{t)dt 


wRW~5RW 


+ 


2  IT 

/“ RW 

w(3ujrw~2Srw) 

2wRW  (w RW  ~SRW  ) 


MzdRW(t)dt 


(67) 


Approved  for  public  release;  distribution  unlimited. 
29  of  34 


S3  to 


Again,  let  C\  =  {ujrw  —  $Rw)t  and  C2  =  {urw  +  &rw)  t  +  ^RW-  Substituting  Equations  10  and  17 
into  Equation  67,  while  using  Equations  10  and  17,  yields 


w 


2vr 


ArW(ujrw  ~  8rwY 


W  /  " RW~SRW  .  9  i, 

—tiL—  /  sin  C\ dt 

2  Jo 


—kpAxp  /  'lM  l>M  sin2  C\  sin  {Arw  cosCi  +  i]rw}  dt 


-kLy^p  /  RW  RW  sin2  C\  cos  {ARW  cos  C\  +  rjRW}  dt 


-kLxYpP  cos  a  RW  RW  sin2[(wi?u/  -  $Rw)t]  sin  {ARW  cos  Ci  +  i]rw}  dt 


7  wp  ■ 

—knxA  srna 


°RW  —°RW  .  2 


+  ARw{vrw  +  aR\v)~ 


sin2[(o;w  -  SRW)t]  sin  {ARW  cos  C,  +  r/RWj  dt 


/n{3ulRW^25Rw) 

2“sw("wjw)  .  2  /~7  , , 
sin  C  2  at 


^RW~SRW 


n{3ulRW-2SRw) 


+  kpAxp  /  Rw{ '  sin2  C2  sinjAiju/ cos  C*2  +  ?7/j,vk}  dt 


WRW~SRW 

/A^itw  ~2Srw) 

2uirw(u1rw-srw)  ■  2 

n 


sin2  C2  cos  {A/^vu  cos  C2  +  ??rvk}  dt 


UJRW~°RW 

^{3ivrw-25rw) 

7  WP  /  2"flw(“SH'-%B')  •  2  /-Y  •  r  /t  ,  1  ,, 

+  COS  a  /  Sill  C-2  Sin  {A/Jiy  COS  C2  +  1JRW  s  dt 


U>RW~°RW 

A 3uirw~2Srw ) 

/  WP  ■  I  2wRwipRW~sRw)  ■  2n  ■  (A  r<  ,  l  u 

+  kpxcp  sm  a  /  sm  G2  sm  {Arvv  cos  G2  +  at 


°RW~0RW 


+  (Aru/  +  AArw)  {urw  +  cfrwY 


W  /  URW 

I _A3^fR _ . 

2u1RW  (ujrw~srw  ) 


Al-  /  ;r  \  sill2  CY2dt 

Z  /  ^{^RW-^RW ) 


(68) 


o  /  ^  RW  O 

+  kDAxR  I  n^Rw_2SA  sin  C2  sin  {(ARW  +  AArw)  cos  C2  +  r]Rw}dt 


-  kLy, 


2wrw{ujrw~5rw) 

rz^- 
WP  /  "BW 

cp  /  ’n{3wRW-2SRw) 
2u1RW  iPRW  -  SRW  ) 

2tt 


sin2  C*2  cos  {(Aru^  +  AArw)  cos  C2  +  Pi?w  }  dt 


+  kLx^P  cos  a  /  "W'  sin2  C2  sin  { (Aw  +  AARW)  cos  C2  +  ??w}  dt 

K  /  -k{2,urW~26rw) 

2ulRwiulRW~SRw) 

+  kDx^P  sin  a  /  ^  sin2  C2  sin  { ( Aw  +  AAnw )  cos  C2  +  r]RW  }  dt 

P  J  n(3“RW~2SRW) 

2"bw("bw"jbw) 


/ 
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Note  that  the  integrals  are  of  the  form  given  in  Equations  25-40;  thus,  Equation  68  can  be  written 
as 


M. 


-B  U>RW 


~RW 


2t r 


Arw&rw  -  SrwY 


-  kL^h  -  kDAxRI2  -  kLyWpI3 


—  kLxYp^  cos  al2  —  kDx^p  sin  al2 


+  A?rw{lorw  +  &rwY 


-  kL^I41 


+  ko AxR  (/,5i  +  75 2)  -  kLUcp  {hi  +  hi) 


+  kLx^p  cos  a  (J51  +  152)  +  kDxYp^  sin  a  (75 1  +  752) 


wp. 


cp 


+  (urw  +  <trw)2  |  -  kL~h 2  +  koAxR  (753  +  154) 


-  kLy^p  (763  +  764)  +  kLx^P  cos  a  (753  +  754)  +  kDx]YpP  sin  a  (753  +  754) 


„wp 


WP , 


(69) 


Substituting  the  results  for  the  definite  integrals,  7i-7g,  and  simplifying  gives 


-jrj-B  ArWujrw  {lorw  -  5rw) 


M, 


zRW 


2 

„WPr 


Ji{Arw)  sin tjrw  {~kD  (sin  ax^p  +  AxR) 


W  Jr  7  VV  Jr  7 

—  cosaxcp  ItL  —  ycp  kl  — 


wpu_  cos  Vrw  \  _  wkLARW 
sin  r/Rw  J  4 


ko Arw^rw  (urw  +  (trw)  /  .  wp  ,  a  b\  (  T  t  a  \  ■ 
- - - (sm  ax"  +  Axb )  I  Ji  (Arw )  sin  rjRW 


+ 

-  77i  (ARW0  cos  r/_Riy 

_  kLARwURw  {urw  +  ctrw ) 

4 

-  Hi{Arw)  costirw}  +  {7i(ARw)  cosrjRw  +  H\{ARW)  smr]Rw  [>  + 


—  cos  ax 


WP 

Cp 


Ji(Arw)  sin  t)rw 


wArw 


.  ko  (Arw  +  AArw)urw  (wrw  +  PRW)  s  .  wp  .  A  B\ 
H  -  (sin  axcp  +  l±xr) 

*  (Ji(Arw  +  AArw)  sin  r/Rw  +  H\(Arw  +  A  Arw)  cos  t/rw) 
kL  {Arw  +  A  Arw)  x>rw  {^rw  +  °rw) 


—  cos  ax 


WP 

Cp 


*  {Ji{Arw  +  A  Arw)  sinrjRw  +  Hi{Arw  +  A  Arw)  cos  tjrw} 

+  UcpP  {Ji{Arw  +  A  Arw)  cosr/Rw  —  Hi(Arw  +  A  Arw)  sin  t/rw} 
w  {Arw  +  A  Arw)  \ 

4  J 


(70) 
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Following  a  similar  procedure  for  the  left  wing,  it  can  be  shown  that 


Mzlw  ~  ~ 


Alw^lW  {^LW  —  $lw) 


Ji{Alw )  sin  ijLW  <  -  kD  (sin  ax^pp  +  Axf) 


WPi  WPi  COS7lLW 

-cosax  kL-y  kL— - 

p  sm  r]LW 

koA^W^LW  (y>LW  +  °lw) 


wkLALW 


(sin  ax[:pP  +  Axf)  I  J\  ( ALW )  sin  rjLW 


-  H\{Alw)  cos i]lw 

kpALW^LW  (w LW  +  &Lw) 


+ 


—  cos  ax . 


WP 

cp 


J\{Alw)  sin r/LW 


-  H\{Alw)  cosiilw  }+  Ucp  {Ji{Alw)  cos  r/LW  +  Hi(Alw)  sin  ijlw}  + 


wALW 


kp  ( Alw  +  A  Alw)  ulw  (ulw  +  °lw)  /  .  wp 

- (sm  ax 


+  Axf) 


4  v - °p 

*  (Ji  ( Alw  +  A  Alw)  sin  t/lw  +  Hi  ( Alw  +  A  Alw)  cos  i]lw) 
kL  ( Alw  +  A ALw)  w lw  (w lw  +  &lw) 


+ 


—  cos  ax 


WP 

cp 


*  {Ji{Alw  +  A  Alw)  sin  t/lw  +  H\(Alw  +  A  Alw)  cos  i]lw} 

+  2/cpP  {Ji{ALw  +  A  Alw)  cos  t/lw  ~  Hi(ALw  +  A  ALw)  sin  r]Lw} 
w  ( ALw  +  AAlw)  \ 

4  > 


(71) 


Note  that  without  split-cycle  frequency  modulation,  there  exists  a  non-zero  cycle-averaged  yawing 
moment  on  each  wing  since  y^pp  ^  0  and  w  ^  0.  When  the  split-cycle  parameters  for  each  wing  are 
zero  and  ujrw  =  ^LWi  the  cycle- averaged  moments  are  opposing  and  balance  one  another.  Since  the 
fundamental  wingbeat  frequency,  ui,  can  be  independently  varied  for  each  wing,  yawing  moments 
can  be  generated  without  varying  the  split-cycle  parameters  or  wing  bias  terms.  Nevertheless,  the 
bias  terms  are  also  present  in  Equations  70  and  71  and  therefore  can  be  used  to  generate  yawing 
moments.  This  is  because  there  is  a  difference  in  the  moment  arm  from  left  to  right. 


VI.  Conclusions 

The  analysis  presented  shows  that  the  use  of  split-cycle  constant-period  frequency  modulation 
with  wing  bias  and  independently  actuated  wings  allows  one  to  manipulate  the  x-body  and  z- 
body  axis  forces.  The  y-body  force  is  also  controllable  using  this  approach,  however,  since  small 
wing  bias  terms  are  needed  to  perform  typical  maneuvers,  the  available  cycle-averaged  y-body  axis 
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forces  are  small.  Rolling,  pitching,  and  yawing  moments  can  also  be  generated  using  the  split-cycle 
technique  with  wing  bias  and  it  appears  that  wing  bias  can  significantly  change  the  cycle- averaged 
pitching  moment.  Given  the  ability  to  manipulate  the  cycle- averaged  body-axis  forces  and  moments, 
untethered  controlled  flight  with  insect-like  maneuverability  appears  to  be  feasible  with  only  two 
physical  actuators,  which  is  supported  by  the  simulation  results  presented  in  Part  II.34 
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